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Abstract. I describe the scenario of molecular cloud (MC) evolution that has emerged 
over the past decade or so. MCs can start out as cold atomic clouds formed by com- 
pressive motions in the warm neutral medium (WNM) of galaxies. Such motions can 
be driven by large-scale instabilities, or by local turbulence. The compressions induce 
a phase transition to the cold neutral medium (CNM) to form growing cold atomic 
clouds, which in their early stages may constitute thin CNM sheets. Several dynam- 
ical instabilities soon destabilize a cloud, rendering it turbulent. For solar neighbor- 
hood conditions, a cloud is coincidentally expected to become molecular, magneti- 
cally supercritical, and gravitationally dominated at roughly the same column density, 
N ~ 1.5 x 10 2 1 cnT 2 w 1OM pc~ 2 . At this point, the cloud begins to contract gravita- 
tionally. However, before its global collapse is completed (~ 10 7 yr later), the nonlinear 
density fluctuations within the cloud, which have shorter local free-fall times, collapse 
first and begin forming stars, a few Myr after the global contraction started. Large- 
scale fluctuations of lower mean densities collapse later, so the formation of massive 
star-forming regions is expected to occur late in the evolution of a large cloud complex, 
while scattered low-mass regions are expected to form earlier. Eventually, the local star 
formation episodes are terminated by stellar feedback, which disperses the local dense 
gas, although more work is necessary to clarify the details and characteristic scales of 
this process. 



1. Introduction 

The molecular gas in the Galaxy exists in giant complexes (giant molecular clouds, 
or GMCs) of masses ~ sizes of several tens of parsecs, and mean densi- 

ties n ~ 100 cm (see, e.g., the review hv lBlitzl[l993h . and contain a large amount 
of substructure (parsec-scale clumps of densities n ~ 10 3 cm -3 , and sub-parsec-scale 
cores of densities n > 10 4 cm -3 ). Molecular clouds (MCs) contain roughly half the 
gaseous mass in the Galaxy, and are the sites of all present-day star formation (SF) in 
the Galaxy. Thus, the study of their origin and evolution is crucial for our understand- 
ing of SF, besides the importance, on its own right, of understanding this fundamental 
component of the ISM. 

The seminal paper by lOortl d 19541) postulated the existence of a cycle (now known 
as the Oort cycle), in which the expanding HII regions around newly formed mas- 
sive stars form shells of cold gas around them. The shells subsequently fragment and 
produce a population of cloudlets, which then grow by coagulation until they become 
gravitationally unstable, at which point they proceed to collapse and form a new gener- 
ati on of stars, starti ng the cycle all over again. In this model, which was later developed 
by iField & Saslaw 0961), clouds were assumed to grow exclusively by coagulation 
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of randomly-moving small clo udlets. The sam e clou d growth process was assumed 
in later models, such as that by iNorman & Silkl dl980h . which differed from the Oort 
model mainly in that the driving agent was considered to be winds from low-mass, 
T-Tauri stars r ather than the ionizing radiation from massive stars, and the model by 
McKee ( 1989) for the SF rate (SFR) regulated by the background photo-ionizing ra- 
diation. The coag ulatio n proce ss implied very long (> 10 8 yr) cloud growth times 
(IScoville & Hershl [l979: Kwan 1979), which were however ruled out on the basis of 
observational evidence b y iBlitz & Shul (Il980l) . These authors proposed instead that 
MCs form and grow by a lParkerl (119661) instability triggered by spiral-arm shocks, and 
have lifetimes ~ 10 7 yr. 

Ever since the times of those early models, MCs have been an o dd component 
of the ISM. Because the y are known to be strongly self-gravitating dLarsonl 1 19811 : 



Mvers & Goodman 



1988) and at significantly higher thermal pressures than the mean 



ISM pressure (e.g.. lMyerslll978), they did not fit in thermal-pressure balance models 
of the ISM, such as that by lMcKee & Ostrikerl (1 19771) . Instead, they have traditionally 
been considered to be in approximate virial equilibrium dLarson[l98ll : lMyers & G oodman 
1988), supported against their self-gravity by either the magnetic field (the so-called 
standard" model o f magnetically regulated SF; see, e.g., the reviews by IShu et al. 



19871 : iMouschoviasI 1 1 99 lb . or by turbulen ce driven by stellar feedback ( t he so-called 
"turb u lent" model of SF; see the reviews by Vazquez-Semadeni et al.l2000l: iMac Low & Klessen 



2004; M cKee & Ostrikerl 120071 : iBallesteros-Paredes et all l2007h . In both cases, the 



gravitational contraction was assumed to be halted by either of the two mechanisms, 
and the clouds were assumed to reach near virial equilibrium. In the last decade or 
so, however, the paradigm about the formation, evolution and structure of MCs has 
changed significantly, and in the remainder of this review I will discuss this emerging 
new view. 



2. Birth and infancy 

2.1. Observational and numerical evidence on the clouds' origin 

As mentioned above, GMCs and their substructures had traditionally bee n thought to be 
in vir ial equilibrium. However, recent observations of GMCs in the LMC dKawamura et al.l 
20091) suggest that the clouds are undergoing an evolutionary process, in which both 



their mass and their SF activity increase in time, going in ~ 25 Myr from masses 
M ~ 10 4 5_5 M o and virtually no massive-SF to M ~ 10 5 - 5_6 M Q an d a population of 
cluster s and HII regions. This is consistent with the conclusion by [Engargiol a et al. 



(2003) that the GMCs in M33 are being assembled rapidly from the atomic component, 
with a prompt onset of SF. 

Similar conclusions had been reached previously from numerical studies. Numer- 
ical simulations of the ISM at the kpc scale w ith turbulence driven by stellar feedback 
(supernova explosions, expanding HII regions ;lBania & Lvonlll98~ol:IChiang & Bregman 



1988; Rosen et al. 1993; Rosen & Bregman 1995; Vazquez-Semadeni et al. 1995; Passot et al. 



1995; Korpi et al. 1999; B allesteros-Paredes et al . 1999a; de Avillez 2000; d e Avillez & Breitschw erdt 



2004j, [2005) showed that compressive motions driven by large-scale gravitaional insta- 



bilities in the diffuse ISM or by the global turbulence are able to form clouds on short 
timescales, essentially given by the turbulent crossing time accross the distance nec- 
essary to collect the material that eventually reaches the cloud. This is facilitated by 
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the presence of cooling, which renders th e medium highly compressible, even when no 
thermal instability (cf. 32.2.11) is present (Vazquez-Se madeni et alJfl996D . Within such 
a dynamic scenario. iBallesteros-Paredes et all dl999bl) remarked that the clouds should 
not be considered as isolated objects, because a significant mass flux is expected to ex- 
ist accross their boundaries, since the clouds are being assembled from material from 
the outside. 



2.2. The physical processes 
2.2.1. Instabilities galore! 

The scenario of cloud assembly by convergent motions in the diffuse ISM (e ither driven 
by turb ulence or by large-scale instabilities) was formulated analytically by Hartman n et al 
d200ll) . who pointed out that the c olumn density of cold atomic hydrogen necessary for 
H2 an d CO molecules to form is (Ivan Dishoeck & Black! 1 19881 : Ivan Dishoeck & Blake 
1 19981) 

A'h ~ 1-2 x 10 21 cm" 2 , (1) 



corresponding to Ay ~ 0.5-1, is very sim ilar to the value nece ssary for the same gas to 
become gravitationally unstable (see also lFranco & Cox1[l9 86), and to the column den- 
sity necessary for the gas to become magnetically supercritical, at a typical interstellar 
magnetic field strength of B ~ 5/uG. This implies that when an initially atomic cloud is 
assembled by a convergent velocity field in the diffuse ISM, it should become molecu- 
lar, self-gravitating, and magnetically supercritical at roughly the same time. We now 
discuss this phenomenology in some more detail. 

A fundamental physical ingre dient aiding the formation of dense atomic and molec- 
ular clouds is thermal instability dFieldl 1 1 965 ). which is a consequence of the various 
radiative heati ng and cooling pro cesses operating on the atomic ISM (for a modern dis- 
cussion, see Wolfi re et alJ fl995). The atomic ISM is subject to the so-called isobaric 
mode of this instability: for densities in the range ~ 1-10 cm -3 (T ~ 5000-500 K), the 
gas loses thermal pressure upon a dynamic compression (for a pedagogical discussion, 
see the review by lVazquez-Semadenill2009l) . implying that after the compression it will 
be underpressured with respect to its surroundings, and will continue to be squeezed by 
them until it exits the thermally unstable range, at which point any further compression 
again increases its thermal pressure, until it eventually reaches pressure equilibrium 
with its surroundings, but at a much higher density and lower te mperature. This i s 
the basis of the well-known two-phase model of the atomic ISM (Field et all 1 19691) . 
originating the warm-diffuse and cold-dense phases of this ISM component, which are 
respectively referred to as the warm and cold neutral media (in turn, respectively, WNM 
and CNM). 

The above discus sion assumes that the gas is i n the thermally unstable range to 
begin with. However, iHennebelle & Perault '(1999) showed that transonic compres- 
sions (i.e., of Mach number M s > 1) in the stable warm phase can nonlinearly trig- 
ger a transition to the cold phase, so that cold clouds can be formed out of the stable 
WNM in the presence of transonic turbulence, wh ich is indeed observed in this medium 
dKulkarni & Heilesll 19871 : iHeiles & Trolandll2003l) . Alternativ ely, large-scale (> 1 kpc) 
instabilities in the diffuse medium , such as gravitat ional (e.g.. lElmegreenlfl 994) . Parker 
( Parked IT966h or magneto- Jeans ( Kim et al. |2002), can provide the driving forces for 
these motions. 
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The gas cooled and compressed by this process to form a cloud is subject to a 
large number of dynamical instabilities. It has been long been known that compressed 
layers formed by the collision of gas streams are nonlinearly unstable, meaning that 
the layers become turbulent when the colliding flows h ave sufficiently large velocities 
dHunter et al.lll986l : IStevens et alJll99Alvishiuadll994h . This process is known as the 
nonlinear thin-shell instability (NTSI). Furthermore, the presence of c ooling lowers 
the required inflow velocities to destabilize the layers dPittard et al.ll2005h . Finally, the 
Kelvin-Helmholz and Rayleigh-Taylor instabilities are also expected to operate during 
the formation o f a cloud. The interplay of all these instabilities has been investigated 
numerically by iHeitsch et all d2005l 120061) . In summary, convergent motions in the 
WNM are expected to produce turbulent CNM clouds. 

2.2.2. Evolution of the mass-to-magnetic flux ratio 

Another crucial ingredient in MC dynamics is the magnetic field, and in particular, 
the mass-to-magnetic flux ratio (M2FR). As it is well known, the magnetic field can 
support a cloud against the latter's self-gravity if the M2FR, or, equivalently, the ratio 
of column density to field strength, exceeds some critic al value which, for a cylindrical 
configuration, is given by (iNakano & Nakamura 1978) 



H/B * (4tt 2 G) _1/2 , (2) 

where X is the mass column density. We denote by p the value of the M2FR normalized 
to the critical value. Thus, a magnetically subcritical cloud has p < 1, and can be 
supported by the magnetic fi eld, while a sup e rcritical one has p > 1, and cannot be 
magnetically supported (e.g., IShu et aL 1 11987k lMouschoviaslll99lh . In the "standard" 



model of SF, most clouds were assumed to be strongly magnetically subcritical, and 
thus globally supported by the field. SF was thought to occur on long timescales and 
involving small fractions of the clouds' mass because, in the dense cores, the process 
known as ambipolar diffusion (AD) allows a redistribution of the magnetic flux, so 
that they can continue to contract quasi-statically, until their M2FR finally becomes 
supercritical, and then the cores collapse. 

It is a very common practice to assume that the M2FR is a conserved quantity as 
long as the flow behaves ideally (i.e., AD is negligible). After all, for an isolated cloud 
of fixed mass, the mass is constant by co nstruction , and the magnetic flux is conserved 
by the flux-freezing condition (see, e.g. IShu|[l992l) . However, the M2FR refers to the 



mass within flux tubes, and in general, field lines do not end at the "edge" of a cloud, but 
rather continue out to arbitrarily long distances. In fact, it is well possible that magnetic 
field lines circle around the whole Galaxy. Now, eq. Q implies that a flux tube is 
supercritical beyond an accumulation length given by dMestell 1 19851 : lHartman n et al.1 



2001) 



so that, in principle, the entire ISM is supercritical, at least near the midplane. However, 
in the case of a forming CNM cloud, what is relevant is the M2FR of the dense gas that 
makes up the cloud, since the cloud is up to lOOx denser than its surroundings, and 
thus it is the main source of the self-gravity that the field has to oppose. Thus, in this 
problem, natural boundaries for up to where to measure the M2FR are provided by the 
bounding surface of the dense gas. 
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Accumulation of material is not opposed by the magnetic field if it occurs along 
field lines, and so this can occur freely in t he ISM. When comp r ession s occur at an 
angle with the field, it has been shown by Hennebelle & Perault (2000) that, up to a 
certain angle that depends on the Mach number of the inflows and the field strength, 
the inflow is reoriented along the field lines, and it behaves as if the compression were 
parallel to the lines. Beyond that angle, the compression behaves essentially as if it 
were perpendicular to the field, and the flow bounces off, not forming any cloud. Thus, 
in what follows, we consider the case of a cloud forming along field lines. Equation © 
can be rewritten as 

N CI = 1.45 x 10 21 f-^y cm" 2 , (4) 



5pG 

implying that the column density for a cloud to become supercritical is very similar to 
that for molecule formation to begin, and for the cloud to become gravitationally unsta- 
ble (cf . eq. [T). Thus, we can infer that, as a cloud forms and grows out of a compression 
in the WNM, it should become molecular, supercritical, and gravitationally unstable at 
roughly the same time. 

Figure Q] shows the evolution of the M2FR in two numeric al simulations of clou d 
formation by Vazquez-Semadeni et al. (2011 in prep.; see also [Banerjee et al. 2009), 



illustrating the growth of the M2FR. In these simulations, two cylindrical streams of 
radius 32 pc and length 1 12 pc each, are set to collide against each other at the center 
of a 256-pc numerical box along the x-direction, so that the cloud is a thin cylindrical 
layer perpendicular to the inflows. The M2FR is measured for each line of sight (LOS) 
perpendicular to the cloud, which is defined as the cylindrical volume of radius 32 pc 
and length 20 pc centered at the numerical box center. For each LOS, the M2FR is 
measured as the ratio S/fiy, where B\\ is the field component parallel to the LOS and 
£ is the mass column density along each LOS. Thus, these measurements constitute 
upper limits to the real M2FR (for a detailed discussion, see Vazquez-Semadeni et al. 
2011). The simulations have mean magnetic field strengths of 3 and 4 pG, and p = 0.91 
and 0.68, respectively. That is, both si mulations are globally subcritica l, implying that 



no subregion of it can be supercritical dVazquez-S emade ni et al.l 12005) as long as the 



flow remains ideal. Values of p greater than those of the whole box indicate that the 
measured values overestimate the actual p. Nevertheless, for a real cloud not limited by 
the box size, the magnetic criticality is expected to eventually become > 1. 

These conclusions are consistent with obs ervational results showin g that the CNM 
clouds are in general magnetically subcritical (Heile s & Trolandl 120051. sec. 7), while 



molecular structures appe ar to be critical or supercritical in general (IBourke et aT1l2001 



Troland & Crutcheri l2008) 



2.3. The early stages: thin CNM sheets 

Th e early stages of cloud formati on were investigated analytically and numerically 



by IVazquez-S emadeni et all ( 20061) . determining the structure and physical conditions 



(density, temperature, and expansion velocity of the phase transition front) in the incip- 
ient cloud as a function of the Mach number of the converging gas streams, before the 
dynamical instabilities have time to grow. Figure |2] shows the structure of the cloud (left 
panel) and the dependence of the physical properties of the cloud on the Mach number 
of the inflowing streams {right panel). 

One important implication of the study by Vazquez-Semadeni et al. (2006) is that, 



in its initial stages, the forming cloud is expected to be a thin CNM sheet, since it 
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Figure 1 . Evolution of the normalized M2FR, p, for two numerical simulations by 
Vazquez-Semadeni et al. (2011, in prep.)- See text for details on the simulations. The 
top panels refer to the simulation with (B) = 3pG, while the bottom panels refer to 
the simulation with {B) = 4pG. The left panels refer to measurements performed on 
the entire "cloud" volume, while the right panels refer to measurements perform only 
on those lines of sight for which the number column density satisfies N > 10 21 cm" 2 . 
The solid lines denote the average of p over all LOSs, while the dotted lines denote 
the 3cr deviations from the mean. The growth of p from very small values to near 
unity during the first 10 Myr of evolution is clearly seen, especially in the plots for 
the whole cloud. 

forms at the essentially two-dimensional interface between the colliding streams. After 
a few Myr of evolution, the predicted thin sheet has column de nsities that agree very 
well with the CNM sheets reported by iHeiles & Trolan d (2003). Thus, it is suggested 
that such thin CNM sheets may constitute the earliest phases of MC evolution. Note, 
however, that a GMC may never form if the mass involved in the streaming flows is not 
high enough to attain the column densities necessary for molecule formation. 



3. Maturity: Gravitational contraction and star formation (child bearing) 
3.1. A distribution of collapse timescales 

As discussed in 92.2.21 as a cloud grows, it should become molecular, supercritical, and 
gravitationally unstable at roughly the same time. This result implies that by the time 
a GMC forms, it should be contracting gravitationally, since it cannot be supported by 
the magnetic field because it is already supercritical, and it is not forming stars yet, so 
no additional turbulence that can support the cloud can be injected into it yet. 
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Figure 2. Le/f panel: Schematic illustration of the structure of a cloud formed 
by the collision of oppositely-directed WNM streams. The streams are assumed to 
collide along the horizontal axis. Only the right half of the system is shown, with 
the left stream replaced by a wall. The inflowing gas approaches from the right, 
and shocks at the interface between regions 1 and 2. The shock heating throws the 
post-shock gas out of thermal balance between heating and cooling, and so the gas 
cools as it flows towards the 2-3 interface. At the latter, the gas undergoes a phase 
transition (condenses) to the CNM, and settles into a dense, cold layer (region 3). 
The phase transition front moves outwards from the collision wall at speed Vf . Right: 
Pressure P, number density n, outwards speed of the phase transition front Vf, and 
column density N after 1 Myr of evoluti on, of the thin layer, all as a fun tion of the 
Mach number of the inflow stream. From lVazquez-Semadeni et al.l d2006l) . 



However, the initial turbulence produced by the convergent gas streams has a cru- 
cial role in the subsequent dev elopment of the cloud and its SF activity. Num erical 



simulations by various groups ( Koyama & Inutsuka 2002; I Audit & Hennebellel 120051 . 
iHeitsch et all 20051 : IVazquez-Semadeni et al.l l2006h have shown that the "cloud s" are 
actually a mixture of warm and cold gas (see also Henne beile & Inutsukall200o1) . This 
is illustrated in Vig.^\{left panel), which shows the granular, fractal density structure ob- 
served in a high-resolution (10000 2 zones) numerical simulation by lHennebelle & Audit 
(12007b . The density field has a wide probability density distribution function PDF (Fig. 
|3j right panel), and thus, when gravity is included, t his density PDF implies the exis- 
tence of a wide distribution of free-fall times (Fig. Heitsch & Hartma nn 2008), with 
the bulk of the mass remaining at relatively low densities . This implies that a small 
fraction of the mass in the cloud will collapse, and thus SF will begin, before the global 
collapse of cloud as a whole is completed, on timescales ~ 10 Myr. 

This kind of evolution, which we refer to as hi erarchical gravitational collaps e 
was first observed numerically in the simulations by Vazquez-Semade ni et al. (2007), 
and is very s imilar t o the n otion of hierarchical gravitational fragmentation introduced 
long ago by iHoyld (Il953l) . except that in that proposal, the density fluctuations were 
assumed to be linear, and so they all had essentially the same free-fall time as the large- 
scale cloud, while here they are nonlinear, with shorter timescales than the cloud, due 
to the initial turbulence induced by the very process of assembly of the cloud. 
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Figure 3. Left: Density and velocity fields in a two-dimensional numerical simu- 
lation of colliding flows in a 20-pc numerical box, with a resolution down to a scale 
of 2 x 10~ 3 pc. Note the extremely granular texture of the density field. Right: 
Volume-weighted probability density distribution (PDF) of the density field for this 
simulation (solid line) and for a similar one with lower (half) the resolution (dotted 
line). From lHennebelle & Auditl (120071) . 
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Figure 4. Mass at a given free-fall time as a funtion of time in three different nu- 
merica l simulations of cloud formation by colliding streams by Heits ch & Hartmanr] 
(2008), with the free-fall time axis binned linearly in the top panels, and logarithmi- 
cally in the bottom ones. It is s een that the bulk of the mass remains at low densities, 
with long free-fall times. From Heitsch & Hartmann (2008). 



3.2. Low- and high-mass star-forming regions 

The distribution of density fluctuations produced by the initial turbulence has another 
implication: a clump mass spectrum is produced. From high-resolution simulations, 
Hennebel le & Auditl (120071) reported a power-law shape of the spectrum at high masses 
of the form dN/dM oc M~P, with ft ~ 1.7, in go od agreement with observational deter- 
minations of the mass spectrum for CO clumps (Motte etal.ll 19981 : iKramer et all 19981 : 
Heithausen et al. It is well known that this form of the mass spectrum implies 
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that most of the mass is in the more massive clumps, even though they are less numer- 
ous. 

Of course, the less massive structures are expect ed to be embedded in the larger- 
mass, lower-density ones in a hierarchical fashion (IVazquez-Semadenil 1 19941) . This 
implies that the low-mass structures should be expected to terminate their collapse and 
form stars earlier than the more massive ones, and so a prediction from the present 
scenario is that a large molecular complex should contain a relatively large number of 
low-mass, somewhat older (by a few Myr) star-forming regions, and a smaller number 
of massive regions, formed later, yet containing most of the mass. The latter is consis- 
tent with t he well known res ult that most stars are formed in massive, cluster-forming 
regions (lLada & Ladal 120031) . Of course, this picture neglects triggering of secondary 

SF by previous generations of stars, which complicates the picture; 

This scenario has been quantified in a numerical simulation by lVazquez-S emaden i et al. 



(2009), who selected a typical example of the early-forming, low-mass regions, and the 
most massive cloud formed by colla pse of the cloud complex at lar ge from a numeri- 
cal simulation of cloud formation by Vazquez-Semad eni et all (l2007h. T his simulation 
used a similar setup to that of the si mulations bvlBaneriee et aU (12 009). discussed in 
92.2.21 but with no magnetic fields. IVazquez-Semadeni et al.l (l2009h showed that in- 
deed the masses, sizes, velocity dispersions, and SFRs of the two regions were respec- 
tively consistent with observations of such types of regions. Moreover, they compared 
the distributions of masses, sizes and densities of the dense cores within the massive 
region, showing that they were very sim ilar to the c orresponding distributions for the 
set of cores in the Cygnus-X region by Mo tte et all ( 2007 ). In summary, this analy- 
sis suggested that the formation of low- and high-mass star-forming regions by this 
mechanism is viable. 



4. Old age: Stellar feedback (popping all over) 

Once SF begins in a cloud, the newly formed stars feedback on it via either low-mass 
star outflows, which inject momentum, or ionizing radiation and supernova explosions 
from massive stars. The effect of this feedback on the parent cloud has been extensively 
studied by numerous groups. For a detailed discussion, see the review by Vazquez- 
Semadeni (201 1, in preparation). Here we focus only on whether the feedback is able 
to bring the parent cloud to a quasi-static equilibrium or whether, instead, it ends up dis- 
persing the cloud. This is an issue about which much effort is currently being devoted, 
and because of that, no conclusive answer is yet available. 

I n recent years, simu lations of multiple jets in parsec-scale clumps dNakamura & Li 



20071 : ICarroll et al1l2009l) have suggested that bipolar outflows are sufficient to drive 



and maintain the turbulence in parsec-scale clumps, and to maintain the latter in a near- 
hydrostatic equilibrium. However, in those works the clumps occupied the entire nu- 
merical box, and therefore the simulations lacked the contracting motions of the rest of 
the MC in which they are embedded. 

A simplified model of the effect of massive-star feedback thro ugh HII-region ex - 
pansion on the global evolution of GMCs was performed by Kiumho lz et al. ((2006). 



These authors computed the time-dependent virial balance (not necessarily equilib- 
rium) of a spherical GMC under the effect of its self-gravity and the energy injection of 
its embedded HII regions. In their simplified model, the SFR self-regulates, and causes 
oscillations of the clouds, which are finally dispersed after a few crossing times. 
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Full numerical simulations of HH-region feedbac k in the context of a fully-evolving 
and contracting GMC have been recently performed by lVazquez-Semadeni et al.1 d2010h . 
who again used the setup described in 32.2.21 now in adaptive mesh refinement (AMR) 
simulations which allowed an effective resolution of ~ 0.03 pc, and that included feed- 
back from a single-mass population of stars. These authors found that the feedback 
affects the immediate surroundings of the recently formed stars, thus reducing the SFR, 
but is not capable of reverting the global contraction of the GMC. Similar results have 
be en obtained i n hi gh-resolution simulations at the clump-scale with outflow feedback 
by IWang et al.1 (l2010h . Those authors have found that the accretion that feeds the most 
massive star is not restricted to the core in which the star is being formed directly, but 
rather can be traced out to the scale of the whole clump containing the core, in spite of 
the fact that the outflows distort this flow, reducing the SFE of the system. 

In summary, the available numerical evidence suggests that the scenario of gravi- 
tationally contracting MCs is maintained even in the presence of stellar feedback. The 
latter may be an important source of energy for driving the turbulence at the clump 
(parsec) scale, but has a harder time halting the gravitation al contraction at the sc ale of 
the whole GMC (tens of parse cs). The calculations by iKrumholz et al.1 (120061) and 
IVazquez-Semadeni et al.1 ( 20101) suggest that low-mass clouds are more readily de- 
stroyed by the feedback than more massive ones. Nevertheless, the numerical exper- 
iments performed so far are relatively scarce, and have used a limited set of initial 
conditions. A more comple te coverage of parameter space, such as that conducted by 
iRosas-Guevara et ail (1201 Oh for the variability of the SFR with the initial conditions, 
and with a more accurate modeling of the stellar feedback, is necessary to better under- 
stand the details of cloud dispersal. 



5. Conclusions 
5.1. Summary 

The scenario of MC formation and evolution (under solar neighborhood conditions) 
discussed in the previous sections can be summarized as follows: 

• The route to the formation of a GMC starts with a large-scale, moderately super- 
sonic converging motion in the WNM, which may be driven either by large-scale 
instabilities, the passage of a spiral-arm or supernova shock, or by intermediate- 
scale generic turbulence in the WNM. 

• The compression nonlinearly triggers a phase transition to the CNM, forming 
a large, though not very dense, cold atomic cloud. The earlie st stages of these 



cloud s may constitute the thin CNM clouds reported recently by lHeiles & Troland 



(2003). At later times, a combination of the Kelvin-Helmholz, Rayleigh-Taylor, 
nonlinear thin-shell and thermal instabilities destabilizes the cloud, rendering it 
turbulent and clumpy. 

If sufficient mass is available in the converging flow that a column density ~ 
1.5 x 10 21 cm -2 is reached in the newly formed cloud, then the cloud begins to be 
dominated by self-gravity (rather than by the confining pressure [thermal+ram] 
of its surroundings) and also reaches a high enough extinction (Ay ~ 1) to al- 
low the formation of CO molecules. Moreover, for a fiducial value of the mean 
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Galactic magnetic field of 5//G, such a cloud should be near to becoming mag- 
netically supercritical as well. So, for solar neighborhood conditions, a growing 
cloud should become molecular, self-gravitating and magnetically supercritical 
at approximately the same time. 

• A cloud that has reached such a column density then begins to contract gravita- 
tionally. However, the clumps produced by the various instabilities, which have 
shorter free-fall times than the bulk of the cloud, culminate their collapses be- 
fore the bulk of the cloud does. Star formation thus begins a few Myr after the 
cloud's global contraction has started, but a few Myr before the bulk of the cloud 
cuminates its own collapse. 

• The collapse of isolated, low-mass clumps produces scattered low-mas star-form- 
ing regions, while the collapse of the bulk of the cloud produces high-mass star- 
forming regions. 

• The termination of a SF episode in a cloud is still not fully understood from 
the avilable numerical simulations. More complete coverage of parameter space 
and with more detailed modeling of the feedback is necessary. In any case, the 
attainment of a nearly hydrostatic equilibrium appears very difficult. Instead, it 
appears that a cloud may begin to be destroyed or dispersed locally, while the 
outer layers may still be falling in, establishing a complicated flow with both 
infall and outflow. 



5.2. Implications: hierarchical gravitational contraction 

The evolutionary scenario for MCs described here strongly suggests the ubiquitous ex- 
istence of generalized gravitational contraction in MCs and their substructure, given 
the approximate simultaneity of the onset of gravitational contraction with the onset 
of molecule formation and the attainement of magnetic supercriticality. This implica- 
tion is in fact consistent with a growing body of observational results. In particular, 
principal component analysis of the contributions to the velocity dispersion in MCs 
and their substructure invariably show a "dipolar" main component, indicating that the 
dominant contribution come s from a large-scale ( i .e., at the scale o f the full structure 
observed) velocity gradient dHeyer & Bruntl 120071 : iBrunt et al.ll2009l) which is consis- 
tent with a global contraction or shear of the cloud, but inconsistent with solid-body 
rotation (M. Heyer, pri v. comm., 2010). Also, studies comparing spe cific regions with 



numerical simulations (Hart mann & Burkertll2007klPeretto et al.ll2007f) have shown that 



those regions are well modeled by gravitation ally contracting structures . Finally, recent 



studie s of massive star-forming regions by dGalva n-Madri d et"al] |2009: Csenge ri et al. 



1201 Ol) have provided evidence of the existence of hierarchical accretion flows at mul- 
tiple scales in the regions. The notion of hierarchical gra vitational fragmentation has 
recently been formulated analytically by Fie ld et all ( 20081) . who have proposed the ex- 



istence of a gravitationally driven cascade from the large to the small scales in MCs, 
analogous to a turbulent cascade, except that the quantity being cascaded is mass rather 

than energy. 

Moreover, recent observations bv lHevereta l. (2009) having much higher angular 



and s pectral resolution and higher dynamic range in column density than earlier studies 
(e.g.. lSolomon et aT1ll987l) . suggest that MCs do not, after all, have a roughly constant 
column density E, but rather span up to two orders of magnitude in this variable, and 
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that the velocity dispersion actually scales with size and column density as 

o-v oc (SL) 1/2 . (5) 

Recently, iBallesteros-Paredes et all d2010h have proposed that this scaling is exactly 
what is expected from a hierarchical gravitational cascade, in which, rather than virial 
equilibrium, the governing relation is simply energy conservation during the contrac- 
tion, so that the gravitational and kinetic energies satisfy 



|£ g l = E k , (6) 

from which eq. (f5]) follows directly. It is noteworthy that, at face value, eq. © seems 
to fit the data better than the virial equilibrium condition \E g \ = 2E k , as shown in 
Fig. 12 although large uncertainties, especially in the mass determinations, prevent any 
firm conclusions. Interestingly, note that free-fall implies larger velocities than virial 
equilibrium, contrary to the standard notion that velocities higher than virial imply lack 
of gravitational binding. 




Figure 5. Mas sive dense cores fro m|Gibson et al. (2009, labeled "G"), and clouds 
and clumps from lHeyer elaTl (I2009L labeled "H") in the cr v /R 1 ^ 2 vs. S plane, where 
cr v is the velocity dispersion, R is the region size, and S is the mass column density. 
The straight lines show the loci of virial eq uilibrium, \E e \ - 2E\, and of e nergy 
conservation under free-fall, |£g| — E^. From Ballesteros-Parede s et al.1 (120101) . 



It is important to remark that the hypothesis that MCs might b e undergoing gener- 
alized gravitational contraction is not new. It was fir st proposed by Goldreich & Kwan 
(|l974h . However, it was soon deemed untenable by IZuckerman & Palmerl (| 19741) who 
argued that, if the clouds were in free-fall, then a simple estimate obtained by divid- 
ing the total molecular mass in the Galaxy by the mean free-fall time would imply a 
SFR roughly two orders of magnitude larger than observed. However, this conundrum 



Molecular Cloud Evolution 



13 



may not pose a serious problem if the star-forming activity of the clouds is terminated 
prematurely by the feedback from the first stars formed. More work to understand the 
details of this process is still needed. 

5.3. Final remarks 

The picture of MC evolution described in this review appears promising as a self- 
consistent scenario that connects the dynamics of the ISM at large with the physical 
properties of MCs and their star-forming properties. However, several features still 
need to be worked out in more detail, such as the precise form in which stellar feed- 
back terminates a SF episode and disperses a cloud locally, and at what scales can this 
process be considered in an averaged sense, if at all. Also, more quantitative statistical 
issues, such as the fraction of the time, and of the stellar production of the clouds, is 
spent under the magnetically subcritical and supercritical regimes, respectively. 

Finally, it is important to remark that the present scenario is likely to only apply 
directly to solar neighborhood-like conditions. In particular, it may need modification 
for regions like the molecular ring of the Galaxy, or to galaxies where the disk is mostly 
molecular, as in those cases the atomic component may be absent, at least locally. In 
those cases, it is possible that the entire molecular ring or disk is the equivalent of 
the isolated GMCs we have discussed here, and that the phase transition occurs at the 
boundaries of these regions, both radially and vertically. More work is needed in order 
to assess this. 
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